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Domain- Wall Fermions at Strong Coupling * 

L. Levkova a and R. Mawhinney a 

a Department of Physics, Columbia University, New York, NY, 10027 

The DWF formulation becomes increasingly problematic at gauge couplings for which a -1 < 2 GeV, where the 
roughness of the gauge field leads to increased explicit chiral symmetry breaking (m rcs ). This problem becomes 
especially severe for sufficiently strong coupling where the underlying 4-dimensional Wilson theory is in the Aoki 
phase. We review our attempts to find a suitable modification of the gauge and/or the fermion action which 
would allow the DWF method to work reliably at stronger coupling. 



1. INTRODUCTION 

Domain-wall fermions have the advantage of 
having an exact chiral symmetry in the limit of 
an infinite 5th dimension on the lattice, even at 
a non-zero lattice spacing. In combination with 
the DBW2 gauge action £Q in quenched simula- 
tions it has been shown that the r7i rcs is about 
two orders of magnitude smaller than its value 
for the Wilson gauge action [2]. This effect has 
been attributed to the properties of DBW2 to 
suppress lattice dislocations, which contribute di- 
rectly to the m res , and to make perturbative cor- 
rections small |2I3) . However, in full QCD calcu- 
lations at coarse lattice spacings, for the combi- 
nation of DWF/DBW2, the m rcs is still not sat- 
isfactory small. The quark-gluon system is in the 
Aoki phase where the DWF formulation suffers 
from the violations of chiral symmetry and local- 
ity. Our goal is to find a suitable modification of 
the gauge action and/or the fermion action which 
would allow us to conduct a thermodynamics cal- 
culation with DWF. 



EFFECTS OF THE ADDITION 
AN ADJOINT TERM TO DBW2 



OF 



Our studies of the plaquette distributions of dy- 
namical DWF/DBW2 configurations 0] showed 
that with the increase of the coupling the tails 
of the distributions extend to more negative pla- 
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quette values, which are associated with the ap- 
pearance of greater number of lattice dislocations 
and increased m TCS . These findings gave us the 
idea to modify the DBW2 action by adding an 
adjoint term to it with the expectation that the 
modified action will suppress more strongly the 
negative tails of the plaquette distributions and 
thus it will get rid of the lattice dislocations. The 
modified DBW2 action has the following form: 



Sdbw2+adj 



JL 



(1 -8ci) 



^2 RoTtP^(x) 



where we choose c/ + 2c a = 1, to keep the nor- 
malization for this three term action the same as 
for the plaquette action in the continuum limit. 
The tail shortening effect on the plaquette dis- 
tributions is an expected result of the decreasing 
of the effective coupling for the fundamental pla- 
quettes in the action, when the coefficient c a is 
negative. 

In the quenched case, we chose to match the 
plaquette distributions of two runs, one with the 
modified DBW2 and the other with the original 
DBW2, and compare their m res and lattice spac- 
ing. Figure ^ shows that although the plaque- 
tte distribution of the modified DBW2 run b has 
a shorter negative tail, its m rcs = 0.0156(2) is 
actually an order of magnitude larger than the 
corresponding value of m ICS = 0.00125(3) for the 
unmodified DBW2 run a. Furthermore run a is 
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out of the Aoki phase and run b according to its 
closed spectral flow gap shown on on Figure [3 is 
in that phase. In conclusion, we attempted to 
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Figure 2. Spectral flow of jsDw for run a and b 
from Figure n 



Plaquette 



Figure 1. Plaquette distributions for quenched 
runs a and b. Run a has /3 = 0.87 and a^ 1 k 1.3 
GeV. Run b has /3 = 0.40, c/ = 2.0 and a" 1 « 
0.93 GeV. Both runs have volume 16 3 x 32 and 
L s = 12. 

"cure" the Aoki phase by changing the action at 
the scale of the lattice spacing by adding an ad- 
joint term and shortening the tails of the plaque- 
tte distributions. The negative tail of the plaque- 
tte distribution although probably contributing 
to m res out of the Aoki phase, when the system 
is in that phase, does not play a significant role 
in determining m ros . The extent of the tail of the 
plaquette distribution into negative values is not 
a primary reason for the onset of the Aoki phase 
and although the DBW2 with the added adjoint 
term does reduce it at a given lattice spacing, that 
does not help to push the system out of the Aoki 
phase. 

3. EFFECTS OF A LARGE RECTAN- 
GLE TERM IN THE GAUGE AC- 
TION 

The rectangle term in the gauge action reduces 
the effective coupling for the plaquette term and 
"smoothes" the gauge field on the lattice by intro- 
ducing some amount of non-locality at the scale 
of two lattice spacings. Since we have concluded 



from the previous section that changing the ac- 
tion at the scale of one plaquette does not influ- 
ence the Aoki phase, which is a more long-range 
global phenomenon, we want to investigate the 
effects of adding more of the less-local rectan- 
gle term. In practice we change the coefficient 
Ci = —1.4069 in the DBW2 action with a larger 
negative value. 

We performed a quenched calculation with a 
DBW2-style action with = 0.53 and ci = -2.3, 
which has a' 1 « 1.0 GeV and m res = 0.0035(1). 
The spectral flow for this run on Figure |31 shows 
a gap and some crossings which we believe means 
that this run is out of the Aoki phase even at that 
coarse a. However to use such a large rectangle 
term for physics calculations could have undesired 
consequences, since the short distance physics is 
distorted by the introduced non-locality. 

4. TWISTED MASS TERM IN THE 
DWF DIRAC OPERATOR 

As another attempt to solve the problem posed 
by the Aoki phase we would like to add an irrel- 
evant to the physics term in the DWF action, 
which would regularize the approach to the Aoki 
phase. We want to experiment with the term 
im T "/^T 3 S XtX 'S s , s i , called a "twisted mass" term. 
As an analogue of M5, the twisted mass term 
provides another parameter, m T , which could be 
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Figure 3. Spectral flow of for a quenched 

run with DBW2-style plaquette plus rectangle ac- 
tion where (3 = 0.53 and c\ — —2.3. Volume 
16 3 x 32 and L s = 12. 



independently tuned to hopefully achieve a bet- 
ter localization of the bound to the domain-walls 
light chiral states, which, if possible, could help 
with the Aoki phase problem. We have investi- 
gated the free case analytically to make sure that 
the twisted mass term addition does not break 
down the most important properties of the DWF 
formulation, namely, the existence of light bound 
to the domain walls states. We tested the ef- 
fects of the twisted mass term on the eigenvalues 
and eigenmodes of D T D^ , by running the Ritz 
algorithm on one chosen lattice from the dynam- 
ical DWF/DBW2 @ = 0.75 run. The graph in 
the middle of Figure 0] shows the dependence of 
the lowest eigenvalue of D T D+ on m T . We see 
that for the range m T < 0.8, the eigenvalue is not 
strongly affected by the value of the twisted mass, 
a behavior very similar to the one observed in the 
free theory (shown on the same graph). This is 
encouraging since it shows that m T , as we hoped 
by its design, is an irrelevant parameter for a cer- 
tain range of values, in the sense that it does not 
act as an additional fermion mass in the prob- 
lem. The dependence on m T of the density in the 
fifth dimension of the eigenvector corresponding 
to the lowest eigenvalue can be studied from the 
bottom graph on Figure Obviously the chiral- 
ity and the localization of this mode changes with 
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Figure 4. Up: Spectral flow of the test lattice. 
Middle: Lowest eigenvalue of D T D+ dependence 
on to t in the interacting and the free field case. 
Bottom: Eigenvector density distribution in the 
5th dimension for different values of m T . 



to 7- and we expect that this is true for the rest of 
the eigenmodes. Whether we can manipulate the 
system by varying m T in a way that globally the 
zero-eigenmode condensation in the Aoki phase is 
suppressed, is still to be determined and will be 
a subject of future investigations 

5. CONCLUSIONS 

The addition of an adjoint term to the DBW2 
action in order to reduce further TO re s at strong 
couplings did not give a significant improvement. 
Increasing the amount of the rectangle term in 
the gauge action gives an improvement in m lea at 
coarse scales, but at the possible expense of the 
distortion of short distance physics. 

Finally, we showed that the twisted mass term 
in the DWF Dirac operator does affect the local- 
ization of the eigenvectors, which we hope is an 
evidence that by tuning m T we can achieve bet- 
ter localization of the eigenstates and influence 
the onset of the Aoki phase. 
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